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Abstract 


We study the structure of vorticity and hydrodynamic helicity fields in peripheral heavy-ion 
collisions using the kinetic Quark-Gluon String Model. The angular momentum conservation within 
this model holds with a good accuracy. We observe the formation of specific toroidal structures of 
vorticity field (vortex sheets). Their existence is mirrored in the polarization of hyperons of the 
percent order. 
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I. INTRODUCTION 


fl 


The local violation |l| of discrete symmetries in strongly interacting QCD matter is now 
under intensive theoretical and experimental investigations. The renowned Chiral Magnetic 
Effect (CME) uses the (C)P-violating (electro)magnetic field emerging in heavy ion collisions 
in order to probe the (C)P-odd effects in QCD matter. 

There is an interesting counterpart of this effect, Chiral Vortical Effect (CVE)^ due 
to coupling to P-odd medium vorticity leading to the induced electromagnetic and all 


conserved-charge currents [sj], in particular the baryonic one. 

Another important P-odd observable is the baryon polarizatiion. The mechanism anal¬ 
ogous to CVE (known as axial vortical effect, see { 4 ] and references therein) leads to in¬ 
duced axial current of strange quarks which may be converted to polarization of A-hyperons 
B fl 0. Another mechanism of this polarization is provided by so-called thermal vorticity 
in the hydrodynamical approach [7[. 

The zeroth component of axial current and correspondent axial charge are related to 
hydrodynamical helicity 

H = J dV(v ■ iu), 

being the projection of velocity v to vorticity w = curlv. This quantity is manifesting 
the recently discovered [5J] and confirmed [g] phenomenon of the separation, i.e. its mirror 
behavior with the same magnitudes but different signs in the half-spaces separated by the 
reaction plane. 

The noncentral heavy ion collisions could naturally generate a rotation (global or local, 
both related to vorticity) with an angular velocity normal to the reaction plane, which is their 
generic qualitative feature. It is naturally to expect that angular momentum conservation 
plays an essential role in the defining the quantitative properties of vortical effects. At the 
same time, it remains to be studied to which extent the particles carrying the main part of 
angular momentum participate in the collisions. 

In the current paper we address these problem by performing the extensive numerical 
simulations. We explore the distribution of angular momentum and find that the role of 
participant nucleons is relatively small, albeit noticeable. We study in some detail the 
structure of vorticity field and compare different approaches to polarization calculation. We 
observed the peculiar toroidal ’’tire-like” structure manifesting themselves in the polarization 
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of hyperons. We will explore different approaches to polarization calculation which will lead 
to qualitatively similar results. 


II. ANGULAR MOMENTUM CONSERVATION IN THE KINETIC MODEL 


The natural source of the P-odd observables in heavy-ion collisions is the pseudovector 
of angular momentum. The question immediately emerges whether it is conserved in the 
course of evolution governed by Quark-Gluon String Model (QGSM) [9]. To check this we 
calculated the angular momentum at various time moments of collision taking into account 
both the contributions of participants and spectators. We consider the Au + Au collisions 
with b = 8 fm at \/i = 5 GeV/u typical for future NICA collider. We observed (see Fig.l) 
that the participants carry about 20% of angular momentum and that the total angular 
momentum of participants and spectators is conserved with a rather good accuracy. 



Figure 1: The time dependence of the total (My) and fireball (Mp) angular momenta in Planck 
constant units and that of hydrodynamical helicity ( H ). 
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One may conclude that the angular momentum is under the good control in QGSM 
model. 

We also studied correlations of angular momentum and the hydrodynamic hclicity and 
observed that these quantities vary in accordance with each other. In order to perform 


this comparison and further ca 
following the earlier suggested 


dilations, velocity, vorticity and helicity were determined 
ij] procedure when the respective quantities were properly 
averaged over events and particles within the three-dimensional cells providing the transition 
from the kinetic to hydrodynamic description. 

Let us pass to the corresponding results. 


III. LARGE-SCALE STRUCTURES OF VORTICITY FIELDS 


We start our studies with the qualitative structure of velocity and vorticity fields. 

The general structure of velocity field follows the ’’little bang” pattern which may be 
quantified by the velocity dependence allowing to extract the ’’little Hubble” constant. We 
calculated the dependence of average cell velocity on the transverse distance p = \Jx 1 + y 2 
and found (see Fig.2) that it is consistent with the linear ’’Hubble” law 


< v/c >= Vq/c + Hp. 


( 1 ) 


The Hubble constant H is measured in the units 10 22 sec = 30 fm/c and changes in the 
range 

H = 0.024 -=- 0.028 (/m/c)" 1 . 

It corresponds to the ’’little Universe” lifetime of about 40 fm/c which is only twice larger 
than the collision time. 

Our key observation is that while velocity field represent the ’’little bang” picture, vor¬ 
ticity field form the relatively thin toroidal ’’tire-like” structures (Fig. 3), which emerge in 
the layer (where velocity field changes rapidly) separating the ’’core” and ’’corona” regions 
loi . 11 ] and form the sort of vortex sheet. 

The interesting property of these structures is that, while emerging due to angular mo¬ 
mentum pseudovector M in the non-central collisions they do not ’’remember” the pro¬ 
duction plane and possess the cylindrical symmetry w.r.t. collisions axis z. This may be 
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Figure 2: The cell velocity dependence on the transverse distance. 

observed (Fig.4) by considering the vortex sheet in the case of the particular direction of M 
along the y axis. 

Such behaviour may resemble cyclones appearing at femtoscopic scale. 

Let us now discuss the observable signatures of these nice structures. 

IV. HYPERON POLARIZATION 

We consider hyperon polarization as the observable related to vorticity and hclicity j^j]. 
We shall concentrate mostly on A hyperon production, which has some advantages: they 
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Figure 3: The vortex sheet. 


are produced in large numbers, their polarization may be easily determined in their weak 
decays, and their spin is carried by strange quark. 

We compare the two rather distinct methods of determining the hyperon polarization. 
The first corresponds to ,ts earlier suggested S and explored fl relat.on to the induced 
axial current while the second one follows the procedure based on the thermal vorticity [7[. 
The first method is based on the calculation of strange axial charge 




Nc_ 

27r 2 


d 3 x n 2 (x) 7 2 eu, dj u f. 

< > N C H 

~ 2 ^ ' 


( 2 ) 


In [5] we used the latter equality exploring the mean-value theorem, while here the spatial 
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Figure 4: The vortex sheet for the particular direction of angular momentum. 

variation of strange chemical potential /i is taken into account. To do so, the description 
of kinetic distribution functions by the correspondent equilibrium equation was performed, 
providing the matching of kinetic and thermodynamical descriptions. As a result, the time 
dependence of the distribution of strange chemical potential takes the form represented at 
Fig. 5. The average polarization can be estimated by dividing Q s 5 (J2]) by the number of 
A’s, assuming that the pseudovector of axial current is proportional to the pseudovector 
of polarization, Ql ~ < Ily ,lab >. Selecting the axial charge related to the particles in the 
def ini te rapidity or transverse momentum interval, the respective dependence of polarization 
may be also obtained. 

As the axial charge should be related to the zeroth component of hyperon polarization 
in laboratory frame fig 16 , the transformation to hyperon rest frame should be performed. 
Taking into account that polarization pseudovector should be directed along y axis (as it 
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Figure 5: The time dependence of strange chemical potential. 


has to be collinear to M pseudovector), one get 


yrA,lab _ 

11 0 — 


n o Py 

Ma ’ 


so that the rest frame polarization can be obtained as 

M a 


< n o > = Ql < 


NxPy 


> . 


( 3 ) 


( 4 ) 


The possible violation of positivity (|IIa| > 1) or even the divergence is due to the fact, that 
hyperons with zero y component of the momentum should not have the zeroth component 




















of polarization and therefore should not contribute to Q s 5 . To avoid this problem one may 
instead attribute the factor p y /M to each hyper on in the denominator of fjTjl . Nevertheless 
the comparison (see Fig. 6) of various approaches shows the similar scale and rapidity 
dependence of polarization. 


w,=1, w2=1 
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Figure 6: The rapidity dependence of polarization in helicity-based approach. 


Another approach to polarization is based on the so-called thermal vorticity |7j] . To 
provide the comparison we calculated (see Fig. 7) the thermal vorticity field and respective 
polarization. While a scale of the polarization in the thermal vorticity-based approach is 
several times larger, its rapidity dependence, surprisingly enough, appear to be similar in 
these rather distinct approaches. 
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Figure 7: The rapidity dependence of the multiplicity and polarization in the thermal vorticity- 
based approach. The magnitudes of polarization at different times are rescaled by the identified 
factors. 

V. CONCLUSIONS AND OUTLOOK 

We investigated vorticity and hydrodynamical hclicity in noncentral heavy-ion collisions 
in the framework of the kinetic Quark-Gluon String Model. We have confirmed our earlier 
observation that the vorticity field is predominantly localized in a relatively thin layer (2 -T 
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3 fm ) on the boundary between the participant and spectator nucleons and observed that it 
is forming the specific toroidal structures, which might be considered as vortex sheets with 
the unexpected cylindrical symmetry. They look as cyclones appearing at femtoscopic scale. 

The vorticity and hclicity fields are manifested in the A hyperons polarization. We per¬ 
formed its detailed calculations including the simulations of the strange chemical potential. 
We found that the polarization magnitude may reach a percent level. The comparison with 
the very different approach exploring the thermal vorticity leads to qualitatively similar 
results, although the polarization scale is several times larger. 
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